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Abstract 

Chlorosomes are likely the largest and most efficient natural light-harvesting photosyn- 
thetic antenna systems. They are composed of large numbers of bacteriochlorophylls orga- 
nized into supramolecular aggregates. We explore the microscopic origin of the fast excitation 
energy transfer in the chlorosome using the recently-resolved structure and atomistic-detail 
simulations. Despite the dynamical disorder effects on the electronic transitions of the bac- 
teriochlorophylls, our simulations show that the exciton delocalizes over the entire aggregate 
in about 200 fs. The memory effects associated to the dynamical disorder assist the exciton 
diffusion through the aggregates and enhance the diffusion coefficients as a factor of two as 
compared to the model without memory. Furthermore, exciton diffusion in the chlorosome is 
found to be highly anisotropic with the preferential transfer towards the baseplate, which is the 
next functional element in the photosynthetic system. 
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Chlorosomes are the largest light-harvesting antenna systems found in nature and contain hun- 
dreds of thousands of bacteriochlorophylls (BChls) enclosed by a lipid monolayer.-"^ They can 
capture light and transfer it in a form of molecular excitations - excitons - to a reaction center 
with high efficiency. This feature allows bacteria to survive at very low light intensities of around 
4 yumol photons m~ 2 s -1 . 5,6 Electron microscopy (EM) 7-8 studies have suggested that the BChl 
c, d, or e molecules are organized into rod-like aggregates - rolls. Light energy absorbed by the 
rolls is transferred via a baseplate 9 to the Fenna-Matthews-Olson (FMO) complex, and finally to 
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the reaction center, where exciton dissociation occurs. Time scales of the excitation energy trans- 
fer (EET) from the rolls to the baseplate have been estimated by time-resolved spectroscopy. — - — 
The EET from BChl c to baseplate BChl a in chlorosome from Chloroflexus aurantiacus occurs 
for instance within 10 ps.— ^ A recent study 15 has suggested that intra-chlorosomal EET takes 
place in the time range of 117-270 fs. It has been difficult to obtain detailed structural informa- 
tion of chlorosomes owing to their large structural disorder. Thus the microscopic origin of the 
fast and efficient EET is still unknown. Recent interest in the EET in the FMO complex,— - — 
red algae,— ^ and plants 24 motivates this study and suggests that a theoretical analysis of exciton 
dynamics in the chlorosomes might provide insights into their function as arguably nature's most 
efficient light-harvesting systems. 

Recent understanding of biosynthesis of the BChl c in Chlorobaculum tepidum (C. tepidum) 25 
has led to a well-characterized C. tepidum three-point mutant. 26 This mutant synthesizes a chloro- 
some that contains >95% BChl d and forms regularly-packed and tubular- shaped elements. 8 Gana- 
pathy et al.,- 1 have proposed an atomistic structure of the rolls from the three-point mutant by using 
solid-state NMR and cryo-EM images: Syn-anti monomer stacks are basic building blocks, and 
they are organized into coaxial cylinders to form the tubular-shaped elements. Recent measurement 
of the linear dichroism spectra of individual chlorosomes 28 has confirmed that the chlorosomes 
from the three-point mutant are less disordered than those from a wild-type, which is compatible 
with the highly- symmetric structure proposed by Ganapathy et al. 21 The atomistic structures of 
the roll are shown in Figure 1. The syn monomer and anti monomer stacks, Figure 1 (e), run 
perpendicular to the roll symmetry axis; rc-BChl pairs form a ring with rc-fold rotational symmetry. 
These rings are stacked in-parallel so that they constitute a helical cylinder (Figure 1(a)). Several 
concentric cylinders are further combined with a spacing of the order of 2 nm (Figure 1(d)). 

In this work, we combine molecular dynamics (MD) simulations, time-dependent density func- 
tional theory (TDDFT) excited-state calculations, and quantum stochastic propagation of the exci- 
ton dynamics to characterize exciton diffusion in the system proposed by Ganapathy et al. 21 We 
find that the supramolecular structure of the rolls enables transfer in the circumference direction, 
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which is relevant for transfer to the baseplate. The environment induces the fluctuations in the elec- 
tronic transitions of the BChl in the roll. The memory effects associated to these fluctuations assist 
the exciton diffusion through the roll and can result in the enhancement of the exciton diffusion as 
a factor of two as compared to the model without memory (Markovian model). 

We briefly summarize the theoretical methods for the exciton dynamics. Details of the calcula- 
tions and underlying assumptions of the theory will be presented elsewhere. Each Bchl d is treated 
as a two-level system of ground and Q y excited states. The time-dependent excitonic Hamiltonian 
coupled to the environmental degrees of freedom is given by 

N N 

H{t) = Y J {(e m ) + Ae m (R(t))}\m)(m\ + £ V mn \m){n\, (1) 

where the \m) denotes a state where excitation is localized at the m-th BChl (site) and all other 
pigments are in the ground state. R(0 is a nuclear configuration at time t. Here, Ae m (R(i)) = 
e m (R(0) - (e»iX e m(R(0) is an excitation energy of m-th site in the nuclear configuration R(0, and 
the bracket () denotes an average over the nuclear configurations. The electronic coupling between 
m-th and n-th sites V mn is calculated within the point-dipole approximation. It is assumed that 
transition dipole moments are oriented along the line connecting A^i to iVni atom in a BChl and 
have 30 Debye 2 as the squared dipole strength estimated from the experimental data. 12 

The time evolution of the exciton system is given by the time-dependent Schrodinger equation: 

i%j^(t))=H{t)\m)- (2) 

Following our previous work on the FMO complex,— we take a stochastic approach to obtain the 
fluctuations of the site energies. Site energy fluctuations are obtained by solving a set of Langevin 
equations with the white noise 

d Ae(0 

-Ae m (t) = —+F m (t), (3) 

at r m 
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where F m (t) is a random force; (F m (t)) = 0, and (F m (t)F m (0)) = 2(Ae%J S(f)/T m . According to 
the eq. 3, the site energy correlation function is (Ae m (t)Ae m (0)) = ^Ae^exp(-t/r m ). Therefore, 
the exciton dynamics governed by eq. 2 is equivalent to the Kubo-Anderson (KA) stochastic 
model.-2L22 Here, a relaxation time of the site energy fluctuation r m is defined as an integration 
of an autocorrelation function (ACF) of the site energy. The density matrix is obtained as an 
ensemble average of these unitary evolutions, p(t) = (1/M) Zf^KOX^KOl- The input parameters 
to this description of the exciton dynamics are (Ae m ), (Aefy, and r m , which are obtained from 
MD/TDDFT calculations. We use a model of ten-stacked rings, Figure 1 (a), for the exciton 
dynamics in the chlorosome, where each ring consists of 60 BChls. Due to the symmetry of the 
roll, we can assume that the statistical properties of the environmental noise are similar for all the 
syn-anti stacked monomers. Therefore, a random sampling of MD/TDDFT trajectories for a single 
dimer of syn and anti BChls can be used for the whole roll. 

To include the contributions from several rolls, we set up a model of three-concentric five- 
stacked rings, Figure 1 (c) and (d), for MD/TDDFT calculations. Rings from different layers 
consist of 60, 80 or 100 BChl molecules, and the total number of the BChl molecules is thus 
1200. A time step of 2 fs was employed, and the temperature was set at 300 K by the Langevin 
thermostat in the MD. The hydrogen atoms were constrained at their ideal positions using the 
SHAKE algorithm. After a 1 ns equilibration run, 20 ps production run was carried out and stored 
every 4 fs. The MD simulations were performed by a CUDA implementation of NAMD. 31 ' 32 
Using the snapshots from the 20 ps trajectory, the Q y excitation energies for the syn and the anti 
monomers were obtained to calculate the averages and correlation functions of the site energies. 
TDDFT excited- state calculations were performed with the long-range corrected hybrid functional 
of Becke (a>B97) 33 and 6-3 1G basis set in the presence of external point charges of all other 
molecules, as implemented in the Q-CHEM quantum chemistry package. 34 

The average site energies were calculated to be 15243.0 and 15319.4 cm -1 for the syn and 
the anti monomers, respectively, and the corresponding standard deviations were 532.2 and 556.5 
cm -1 . Although the syn and the anti monomers are not in same environment, the difference in the 



5 



averages and standard deviations of the site energy are small. An absorption maximum was esti- 
mated by diagonalizing the Hamiltonian without Ae m (t) and was obtained as 13220.6 cm" 1 . This 
value is in reasonable agreement with the experimental absorption maximum of 13888.9 cm" 1 .— 
The ACFs and spectral densities for the syn and the anti monomers are shown in Figure 2. The 
spectral densities were calculated as a weighted cosine transform of the correlation function 



The calculated relaxation times for the syn and the anti monomers are 8.80 and 4.51 fs, respec- 
tively. It is remained to be elucidated whether this difference has some physical importance or it 
stems from the approximations of the model. The spectral densities show the vibrational modes 
in the roll coupled to the excitation energies. Strong peaks around 1600 to 2000 cm -1 can be at- 
tributed to internal bond stretching in the porphyrin macro cycle of the BChl. Characteristic peaks 
are found in the low-frequency region, which are ascribed to intermolecular interactions arising 
from the closed-packed structure. Reorganization energies were calculated by Er = du)J(a>)/a> 
and are 266.5 and 298.4 cm -1 for syn and anti monomers, respectively. 

We then run exciton dynamics simulations by solving eqs. 2 and 3. Initial conditions of the 
wavefunctions were chosen as a localized state on an anti monomer (Figure 1 (b)), and Ae ra (0) were 
sampled from the Gaussian distribution. The density matrix was calculated by averaging 1000 re- 
alizations of the stochastic processes. Figure 3 shows population dynamics obtained by the exciton 
simulations. Excitons delocalize over the entire roll at the time scale of approximately 200 fs. 
The population spreads faster in the circumference than in the coaxial directions. This anisotropic 
diffusion is due to the difference in the electronic coupling constants. Table 1 shows the electronic 
coupling constants between the initially excited monomer and its neighboring monomers. The 
electronic coupling between syn-anti monomers are larger than those between monomers from 
neighboring rings. 

In order to quantify the exciton diffusion, second moments of the exciton propagation in the 
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circumference and the coaxial directions (see Figure 1 (a)) were calculated and are shown in Figure 
4 (a). The second moments initially scale quadratically and then linearly in time. The transition 
from the initial ballistic regime to the diffusive regime is observed at around 20 fs, which is simi- 
lar to high exciton mobility J-aggregates. 37 The second moment in the circumference direction is 
about four times larger than in the coaxial direction. This anisotropic diffusion may play a role 
in the efficient EET: as it is proposed that the symmetry axis lies parallel to the baseplate, 9 the 
faster diffusion in the transverse direction implies the efficient EET to the baseplate. Furthermore, 
the present method is compared with the Haken-Strobl-Reineker (HSR) 38 ' 39 and the classical hop- 
ping 40 models. In the HSR and the hopping models, a pure dephasing rate for m-th site y m was 
obtained as y m = 2(Ae^T m /Ti.— The KA model gives the second moment twice as large as those 
by the HSR, leading to memory-assisted diffusion (MAD). The difference between the KA and 
the HSR is in the memory of the bath fluctuations (KA: (Ae OT (f)Ae m (0)> oc exp(— t/r m ) vs HSR: 
(Ae m (0Ae m (0)) oc 6(t)). This comparison indicates that the memory of the bath fluctuations en- 
hances the diffusion. The diffusion of the excitation energy is directed by the supramolecular 
arrangement and enhanced the memory of the bath fluctuations, which, according to the results in 
this paper, might be very essential for the fast and efficient EET. These results are complementary 
to our findings that non-Markovianity is also near maximal for the FMO complex.— 

Finally, we investigate the dependence of the exciton diffusion on static disorder. Origins of 
the static disorder are, for example, inner structure disorder inside the chlorosome and existence 
of different homologues of the BChls. Here, the static disorder was incorporated by introducing 
Gaussian random shifts in the average site energies. Diffusion coefficients in the circumference 
direction were calculated as a function of a standard deviation of variable static disorder parameters 
and are shown in Figure 4 (b). The HSR model gives the same diffusion constants as the hopping 
model. Since rate constants in the hopping model were calculated with (Ae ra (f)Ae m (0)) oc 6(t),— 
their difference is only in the initial ballistic propagation up to 20 fs. The KA models gives about 
twice as large diffusion coefficients as those from the HSR up to the static disorder of 400 cm -1 , 
while the memory effects become small in the presence of the large static disorder. In general, 
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the static disorder broadens an absorption band, and it will be favorable for the chlorosomes to 
capture a wider range of light energy. On the other hand, presence of the static disorder leads to 
the slower diffusions. The MAD may be a mechanism to overcome negative effects of the static 
disorder inherent in the chlorosome, whiles preserving the positive effects of static disorder, both 
working in concert towards the most productive EET. 

We have presented the theoretical study on the exciton dynamics in the chlorosome. We have 
shown that the exciton diffusion is highly anisotropic with the faster direction toward the baseplate. 
The memory effects of the bath fluctuations increase the exciton diffusion coefficients by two times 
for broad range of the static disorder. We have proposed that the supramolecular arrangement and 
the memory-assisted diffusion are the microscopic origin of the fast and efficient excitation energy 
transfer in the chlorosome. The present methods can be applied to more complex systems, such as 
a system that includes many rolls and/or a baseplate. Work in such a direction will provide a further 
microscopic insight into the chlorosomes as well as a design principle of artificial light-harvesting 
systems. 42 

Acknowledgement 

The authors would like to thank Prof. Huub J. M. de Groot for fruitful discussions and the donation 
of the syn-anti chlorosome configuration template. We further appreciate Stephanie Valleau for 
very useful discussions. T. F. and A. A.-G. acknowledge support from the Center for Excitonics, 
an Energy Frontier Research Center funded by the US Department of Energy, Office of Science and 
Office of Basic Energy Sciences under award de-sc0001088 J. C. B. acknowledges support from 
Welcome Trust UK. S. K. S. and A. A.-G. also acknowledge Defense Threat Reduction Agency 
grant HDTRA1- 10- 1-0046. Further, A. A.-G. is grateful for the support from Defense Advanced 
Research Projects Agency grant N66001-10-1-4063, Camille and Henry Dreyfus Foundation, and 
Alfred P. Sloan Foundation 



8 



References 

(1) Olson, J. M. Photochem. and Photobiol. 1998, 67, 61-75. 

(2) Blankenship, R. E.; Olson, J. M.; Miller, M. In Anoxy genie Photo synthetic Bacteria; Blanken- 
ship, R. E., Madigan, M, Bauer, C. E., Eds.; Springer Netherlands, 2004; pp 399-435. 

(3) Oostergetel, G. T.; van Amerongen, H.; Boekema, E. J. Photosynth. Res. 2010, 104, 245-255. 

(4) Hohmann-Marriott, M. E; Blankenship, R. E. Annu. Rev. Plant Biol. 2011, 62, 515-548. 

(5) Overmann, J.; Cypionka, H.; Pfennig, N. Limnol. Oceanog. 1992, 37, 150-155. 

(6) Beatty, J. T.; Overmann, J.; Lince, M. T.; Manske, A. K.; Lang, A. S.; Blankenship, R. E.; 
Dover, C. L. V.; Martinson, T. A.; Plumley, F. G. Proc. Natl. Acad. Sci. U.S.A. 2005, 102, 
9306-9310. 

(7) Staehelin, L. A.; Golecki, J. R.; Fuller, R. C.; Drews, G. Arch. Microbiol. 1978, 119, 269-277. 

(8) Oostergetel, G. T.; Reus, M.; Chew, A. G. M.; Bryant, D. A.; Boekema, E. J.; 
Holzwarth, A. R. FEBS Lett. 2007, 581, 5435-5439. 

(9) Pedersen, M. 0.; Linnanto, J.; Frigaard, N.-U.; Nielsen, N. C.; Miller, M. Photosynth. Res. 
2010, 104, 233-243. 

(10) Savikhin, S.; Zhu, Y; Blankenship, R. E.; Sturve, W. S. J. Phys. Chem. 1996, 100, 3320- 
3322. 

(11) Psencik, J.; Polfvka, T.; Nemec, P.; Dian, J.; Kudrna, J.; Maly, P.; Hala, J. J. Phys. Chem. A 
1998, 102, 4392-4398. 

(12) Prokhorenko, V. I.; D. B. Steensgaard, D. B.; Holzwarth, A. R. Biophys. J. 2000, 79, 2105- 
2120. 

(13) Psencik, J.; Ma, Y. Z.; Arellano, J. B.; Hala, J.; Gillbro, T. Biophys. J. 2003, 84, 1161-1 179. 

9 



(14) Martiskainen, J.; Linnanto, J.; Kananavicius, R.; Lehtovuori, V.; Korppi-Tommola, J. Chem. 
Phys. Lett. 2009, 477, 216 - 220. 

(15) Martiskainen, J.; Linnanto, J.; Aumanen, V.; Myllyperkio, P.; Korppi-Tommola, J. Pho- 
tochem. and Photobiol. 2012, 88, 675-683. 

(16) Engel, G. S.; Tessa, R.; Elizabeth, L.; Ahn, T. T.-K.; Mancal, T; Cheng, Y.-C; Blanken- 
ship, R.; Fleming, G. Nature 2007, 446, 782-786. 

(17) Panitchayangkoon, G.; Hayes, D.; Fransted, K. A.; Caram, J. R.; Harel, E.; Wen, J.; Blanken- 
ship, R. E.; Engel, G. S. Proc. Natl. Acad. Sci. U.S.A. 2010, 107, 12766-12770. 

(18) Mohseni, M.; Rebentrost, P.; Lloyd, S.; Aspuru-Guzik, A. J. Chem. Phys. 2008, 129, 174106. 

(19) Rebentrost, P.; Chakraborty, R.; Aspuru-Guzik, A. J. Chem. Phys. 2009, 131, 184102. 

(20) Rebentrost, P.; Aspuru-Guzik, A. J. Chem. Phys. 2011, 134, 101103. 

(21) Shim, S.; Rebentrost, P.; Valleau, S.; Aspuru-Guzik, A. Biophys. J. 2012, 102, 649-660. 

(22) Collini, E.; Wong, C. Y; Wilk, K. E.; Curmi, P. M. G.; Brumer, P.; Scholes, G. D. Nature 
2010, 463, 644-647. 

(23) Hossein-Nejad, H.; Curutchet, C.; Kubica, A.; Scholes, G. D. J. Phys. Chem. B 2011, 115, 
5243-5253. 

(24) Schlau-Cohen, G. S.; Ishizaki, A.; Calhoun, T. R.; Ginsberg, N. S.; Ballottari, M.; Bassi, R.; 
Fleming, G. R. Nat. Chem. 2012, 4, 389-395. 

(25) Chew, A. G. M.; Bryant, D. A. Annu. Rev. Microbiol. 2007, 61, 1 13-129. 

(26) Chew, A. G. M.; Frigaard, N. U.; Bryant, D. A. J. Bacteriol. 2007, 189, 6176-6184. 

(27) Ganapathy, S.; Oostergetel, G. T; Wawrzyniak, P. K.; Reus, M.; Chew, A. G. M.; Buda, F; 
Boekema, E. J.; Bryant, D. A.; Holzwarth, A. R.; de Groot, H. J. M. Proc. Natl. Acad. Sci. 
U.S.A. 2009, 106, 8525-8530. 

10 



(28) Furumaki, S.; Vacha, R; Habuchi, S.; Tsukatani, Y.; Bryant, D. A.; Vacha, M. J. Am. Chem. 
Soc. 2011,733, 6703-6710. 

(29) Anderson, P. W. J. Phys. Soc. Jpn. 1954, 9, 316-339. 

(30) Kubo, R. J. Phys. Soc. Jpn. 1954, 9, 935-944. 

(31) Phillips, J. C; Braun, R.; Wang, W.; Gumbart, J.; Tajkhorshid, E.; Villa, E.; Chipot, C; 
Skeel, R. D.; Kale, L.; Schulten, K. J. Comp. Chem. 2005, 26, 1781-1802. 

(32) Stone, J. E.; Phillips, J. C; Freddolino, P. L.; Hardy, D. J.; Trabuco, L. G.; Schulten, K. J. 
Comp. Chem. 2007, 28, 2618-2640. 

(33) Chai, J.-D.; Head-Gordon, M. J. Chem. Phys. 2008, 128, 084106. 

(34) Shao, Y. et al. Phys. Chem. Chem. Phys. 2006, 8, 3172-3191. 

(35) Damjanovfc, A.; Kosztin, I.; Kleinekathofer, U.; Schulten, K. Phys. Rev. E 2002, 65, 031919. 

(36) Olbrich, C.; Striimpfer, J.; Schulten, K.; Kleinekathofer, U. J. Phys. Chem. Lett. 2011, 2, 
1771-1776. 

(37) Valleau, S.; Sakin, S. K.; Yung, M.-H.; Aspuru-Guzik, A. J. Chem. Phys., in press. 

(38) Haken, H.; Reineker, P. Z Phys. 1973, 249, 253-268. 

(39) Haken, H.; Strobl, G. Z. Phys. 1973, 262, 135-148. 

(40) Ern, V.; Suna, A.; Tomkiewicz, Y; Avakian, P.; Groff, R. P. Phys. Rev. B 1972, 5, 3222-3234. 

(41) Breuer, H.-P; Petruccione, F. The Theory of Open Quantum Systems; Oxford University 
Press: New York, 2002. 

(42) Scholes, G. D. S.; Fleming, G. R.; Olaya-Castro, A.; van Grondelle, R. Nat. Chem. 2011, 3, 
763-774. 

11 



Table 1: The electronic coupling constants between the initially excited monomer (1) and 
its neighboring monomers (2-7) in the ten-stacked rings. An index in the row refers to the 
monomer labeled by the same index in Figure 1 (b) 
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Figure 1: (a) A model of ten-stacked rings for the exciton dynamics, (b) Initially excited monomer 
(red) and its neighboring monomers in the ten-stacked rings, (c) Side view and (d) top view of a 
model of three-concentric five-stacked rings for the MD/TDDFT calculations, (e) Syn (red) and 
anti (blue) monomers in the three-concentric five-stacked rings, the site energies of which were 
obtained from the TDDFT excited- state calculations. 
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Figure 2: (a) Autocorrelation functions (ACF) of the site energies and (b) spectral densities for the 
syn (red) and the anti (blue) monomers. Low-frequency regions of the spectral densities are shown 
in the inset. 
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Figure 3: Population dynamics at selected times. Pink spheres represent Mg atoms of the BChl 
molecules, and the size of the Mg atom specifies the population value of the BChl molecule. 
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Figure 4: (a) Second moments obtained from the Kubo-Anderson (black), the Haken-Strobl- 
Reineker (HSR) (red), and the classical hopping (blue) models. Solid and dotted lines refer to 
the second moments in the circumference (Mll) and coaxial (Mzz) directions (see Figure 1(a)), 
respectively, (b) Exciton diffusion coefficients in the circumference direction as a function of a 
standard deviation of the static disorder. 
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